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A careful analysis of the gravitational geon solution found by Brill and 



Q\ ■ Bartle is made. The gravitational wave expansion they used is shown to be 

consistent and to result in a gauge invariant wave equation. It also results 
in a gauge invariant effective stress-energy tensor for the gravitational waves 
provided that a generalized definition of a gauge transformation is used. To 
leading order this gauge transformation is the same as the usual one for grav- 
itational waves. It is shown that the geon solution is a self-consistent solution 
to Einstein's equations and that, to leading order, the equations describing the 
geometry of the gravitational geon are identical to those derived by Wheeler 
for the electromagnetic geon. An appendix provides an existence proof for 
geon solutions to these equations. 



1 



Typeset using REVTgX 



I. INTRODUCTION 



Brill and Hartle HI] , BH, developed a very useful method for finding approximate solutions 
to Einstein's equations that correspond to high frequency gravitational waves propagating 
in a background geometry, which is created by the average stress-energy of the waves them- 
selves. In their paper they applied this method to the case of a static spherically symmetric 
background geometry and found that gravitational waves can remain confined in a region 
for a time much longer than the region's light-crossing time. This so-called gravitational 
geon is generated by a large number of high frequency, small amplitude gravitational waves. 
The time average of the curvature due to these waves creates the background geometry of 
the geon, and this background geometry traps the waves for a long time in a region of space 
called the "active" region. [] The BH solution is important because it serves as an example in 
which the gravitational field both creates and responds to its own effective stress-energy. It 
is also an example of a nontrivial (approximate) solution to the vacuum Einstein equations 
that has no curvature singularities. 

There are three reasons for analyzing the solution found by BH in more detail than was 
done in their original paper. The first is that BH did not investigate the question, is the 
stress-energy tensor they used conserved or gauge invariant? To have a self-consistent set of 
equations it is necessary that the effective stress-energy tensor have these properties. Second, 
the validity of the geon solution found by BH has been questioned by Cooperstock, Faraoni, 
and Perry ||||. They claim that the thin shell approximation used by BH to describe the 
active region is invalid, and that gravitational geons due to high frequency, large angular 
momentum waves do not exist. The third reason for analyzing the BH solution in more 
detail is that, because they used a thin shell approximation, BH did not determine the form 
of the geometry within the active region. The active region is the only region where the 
gravitational waves have a significant size, and it is the region where spacetime is most 
strongly curved. It is clearly important to know the geometry of the active region, if one 
wishes to learn anything about the details of the geon solution. 

In this paper a thorough analysis of the geon solution found by BH is presented. A self- 
consistent expansion of the metric and curvature tensors is found. Using results recently 
obtained regarding effective stress-energy tensors for gravitational waves ||, it is argued 
that, to leading order, the wave equation and stress-energy tensor are gauge invariant. It 
is also shown that the stress-energy tensor is conserved to leading order with respect to 
the background geometry. The wave and backreaction equations are explicitly derived in 
the high frequency and large angular momentum limits. It is shown that in and near the 
active region these equations can be cast in a form which is mathematically identical to 
the equations derived by Wheeler for the electromagnetic geon. Thus, to leading order, 
the background geometry of the gravitational geon found by BH is identical to that of the 



1 Gravitational geons are analogous to the original electromagnetic geons of Wheeler Q, which are 
virtual gravitationally bound states of electromagnetic energy. As such, geons have a finite lifetime 
and are not true non-radiative solutions. Gibbons and Stewart Q have shown that non-radiative 
geons, or other exactly periodic solutions, cannot exist in Einstein's theory. 
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electromagnetic geon found by Wheeler .0 This means that the details of the active region as 
discussed by Wheeler are identical (mutatis mutandis) to those of the gravitational geon 
found by BH. There is of course a difference in the wave equations. However, the effective 
time averaged stress-energy of the gravitational waves is, to leading order, identical to that 
of the electromagnetic waves. 

In Sec. II a review of the BH solution is given. In Sec. Ill a self-consistent expansion 
of the wave equation and the effective stress-energy tensor for the gravitational waves is 
presented. Explicit leading order equations are derived in Sec. IV, and it is shown that 
in the active region the equations to leading order are those found by Wheeler for the 
electromagnetic geon. Extension of the calculation beyond the leading order is also discussed. 
Our conclusions are summarized in Sec. V. 

II. THE BRILL-HARTLE SOLUTION 

The gravitational geon found by BH is a solution to Einstein's equations that consists 
of gravitational waves propagating on a static spherically symmetric background, which is 
created by the waves. There is a large number of waves, each with a small amplitude, a high 
frequency, and a large angular momentum. The waves have different angular orientations 
and somewhat different frequencies. The stress-energy of the waves is significant in and 
near a spherical shell called the active region, and insignificant elsewhere. The solution is 
self-consistent in that the geometry produced by the stress-energy of the waves traps them 
in the active region for a long time. 

We begin by reviewing the BH gravitational wave expansion and the BH definition of 
an effective stress-energy tensor for the gravitational waves. Consider a separation of the 
metric into a background part, 7^, and a perturbation, h^, 

9nv = lnv + h^ u . (1) 

The Einstein tensor can also be divided into a part describing the curvature due to the 
background geometry and that due to the perturbation, 

G^igap) = G^^ap) + AG F (7 a/3 , h a p) . (2) 

It is important to note that AG^ is defined by this equation. 

One way to specify the separation (1) is to use some smoothing or averaging procedure 
< > acting on G^ u , and demand that this procedure not affect the value of G^ u {j a f3). In 
addition, one demands both the exact and the averaged Einstein equations for vacuum. This 
leads to 

AG fa ,(>y,h) = <AG fU ,( r ,,h)> (3a) 
G^) = -<AG lw (j,h)> . (3b) 



2 It is already clear from their paper that the geometry BH found in the regions exterior to the 
active region is exactly the same as that found by Ernst for the electromagnetic geon. 
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(Here and hereafter the indices of the arguments of G^ v and AG Mi , are suppressed for 
notational simplicity.) The first equation is regarded as a wave equation, which gives the 
behavior of the perturbed geometry. The second equation is the backreaction equation, 
which describes how the background geometry is affected by the perturbed geometry. From 
the backreaction equation it is natural to define an effective stress-energy tensor for the 
gravitational waves as 

< >= — -— < AG F (7, h) > . (4) 

For any valid perturbation expansion AG M!/ can be written in the form 

AG F = A 1 G F + A 2 G F + ... . (5) 

BH obtained an expansion of the form (5) by first expanding AG in powers of h and its 
derivatives. They then considered the high frequency, large angular momentum limit of 
this expansion. Thus their leading order term AiG consists of the high frequency large 
angular momentum limit of the terms in the original expansion which are linear in h. Their 
second order term A 2 G consists of the appropriate high frequency, large angular momentum 
limits of the terms in the original expansion which are quadratic in h.f\ BH used a time 
averaging for their stress-energy tensor. This plus the high frequency of the waves resulted 
in a stress-energy tensor consisting of the time average of the high frequency large angular 
momentum limits of the quadratic terms in the original expansion of AG in powers of h and 
its derivatives. To completely fix the value of the stress-energy tensor they implicitly made 
the choiceQ < AiG^ >= 0. The resulting wave and backreaction equations are 

A 1 G^( 1 ,h) = (6a) 
G>,(7) = -< A 2 G^( 7 ,/i)> . (6b) 

The background metric for the geon solution is static and spherically symmetric, so it 
can be written in the form 

=diag(-e !y , e A , r 2 , r 2 sin 2 6) . (7) 

BH used a variational approach to find their solution, the essence of which is to show that 
the effective averaged stress-energy tensor has vanishing trace to leading order. This suffices 
to define the metric outside the active region, which is near r = a. The value of a depends on 
the ratio of the angular momentum to the frequency of the waves. In the infinite frequency 



It also contains certain terms from the original expansion that are linear in h. However these 
vanish when a time average is taken. 

4 It is shown in Ref.( ||) that Eq. (4) does not completely specify the definition of the stress- 
energy tensor. An extra condition must be imposed to determine uniquely the value of < AG >. 
Although not presented as such in their paper, the condition < AiG >= imposed by BH served 
to fix the value of < AG > in the gauge they worked in. See Sec. III. A for further discussion. 
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and angular momentum limits the active region shrinks to an infinitely thin shell at r = a. 
BH's self-consistent treatment found in this limit that 



e" = ~ e" A = l r<a, 

e 1 ' = e~ A = 1 r > a , 8 

r 

and that the mass of the geon is M = (4/9)a. Examination of Eq. (8) suggests that, in the 
active region, dX/dr is large while du/dr is of order unity.f] 

III. CONSISTENCY OF THE GRAVITATIONAL WAVE EXPANSION 

Several consistency questions must be answered to validate the BH solution to Eqs. 
(6a) and (6b). First, the effective stress-energy tensor of the gravitational waves, Eq. (4), 
must be conserved with respect to the background geometry 7 in order to be a source of this 
background geometry. Second, the solution ansatz for the wave equation must be sufficiently 
general so that this set of equations can be simultaneously solved. Third it is necessary to 
investigate gauge invariance. For the wave and backreaction equations to be consistent, the 
effective stress-energy tensor must be at least approximately gauge invariant if the wave 
equation is. Also BH used a gauge transformation to reduce the number of components of 
h^. Finally it is necessary to show that the expansion of AG used by BH is a valid one for 
the geon solution. 

It is appropriate to note that we are here constructing one particular type of geon, 
as described in the beginning of Sec. II. No attempt is made to discuss other possible 
gravitational geons, such as thick-shell ones 0. To show existence of geon solutions it is 
only necessary to establish that one particular, carefully selected set of waves does form a 
geon. 



A. Conservation and gauge invariance 

Conservation of the effective stress-energy tensor for gravitational waves can be estab- 
lished in a straight-forward manner. If we have a solution of Eq. (3b), then the Bianchi 
identity satisfied by G^ u (-y) implies that the exact stress-energy tensor (4) is conserved with 
respect to the background geometry 7^. Thus the approximate stress-energy tensor is 
similarly conserved to leading order so long as a self-consistent expansion is used for the 
gravitational waves. 

The question of whether the wave equation and the effective stress-energy tensor are 
gauge invariant is much more difficult to answer. It is useful in this regard to discuss some 
general results, first for the wave equation and then for the stress-energy tensor. 



5 We use the term "unity" to mean the power a n of the geon radius a, with n adjusted so "unity" 
has the same dimension as the quantity with which it is being compared. 
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If the background geometry 7 is a solution to the vacuum Einstein equations then it can 
be shown that if an expansion of the form (5) is used with A n G consisting of terms of nth 
order in h and its derivatives, then the wave equation (6a) is exactly invariant under gauge 
transformations of the form 

hfo/ix) h^ v {x) T/Ja(^')^ t i/ 1av{pt)£, ^ l/iv^^, • (9) 

Here £ results from an arbitrary coordinate transformation of the form 

x" = x li + ^ . (10) 

If the same type of expansion is used but the background geometry is a solution of Eq. 
(6b) and if the transformations are restricted so that £ and its derivatives are no larger 
in magnitude than h and its derivatives, then it can be shown that the wave equation is 
invariant to 0(h 2 ) under the above type of gauge transformation. For high frequency, large 
momentum gravitational waves and the same restrictions on £, Isaacson |I| showed that 



the leading order wave equation is invariant under gauge transformations of the form (9) to 
second order in his expansion. It is not hard, using an argument similar to that used by 
Isaacson, to show that the wave equation for the gravitational geon solution is also invariant 
to second order under gauge transformations of the form (9) so long as £ and its derivatives 
are no larger in magnitude than h and its derivatives. Here and hereafter by rath order we 
mean of the same order as A n G(7, h) in a self-consistent expansion of the form (5). 

Gauge invariance of the effective stress-energy tensor is a different matter. For Eqs. 
(6a) and (6b) to be consistent, the stress-energy tensor must at least be gauge invariant to 
second order. However, Isaacson showed that A 2 G is not invariant under transformations 
of the form (9). The only case he found in which the stress-energy tensor and thus the 
backreaction equation (6b) are gauge invariant to second order is the case of high frequency, 
large momentum waves when the averaging is over a region of spacetime which is large 
compared to the wavelengths of the waves but small compared to the scale on which the 
background geometry varies. 

One might hope that in the geon case, where high frequency large angular momentum 
waves are used, that the stress-energy tensor would similarly be gauge invariant to second 
order. However, an explicit calculation using the gauge transformation (17) below shows 
that this is not the case. In fact, even when the background geometry is a solution to the 
vacuum Einstein equations (including the case of the flat space solution), the stress-energy 
tensor is not gauge invariant if time averaging rather than spacetime averaging is used. Thus 
it appears that Eqs. (6a) and (6b), which were implicitly solved by BH to obtain the geon 
solution, are inconsistent when time averaging is used. 

The resolution to this very serious problem is given in Ref. |J. It is as follows: First 
Eqs. (6a) and (6b) must be replaced by the equations that result from substituting Eq. (5) 
into Eqs. (3a) and (3b). The result to nth order is 

AiG + ... + A n G = < AiG + ... + A n G > (11a) 
G{i) = - < A 1 G+... + A n G> . (lib) 

Then a generalized gauge transformation is used. It is arrived at by using the coordinate 
transformation (10) and not allowing the functional form of the background geometry to 
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change under this coordinate transformation, that is g^ v = 7 M! , + h^ u . Then h^ v is given 
implicitly by the equation 

7^(x) + hp,(x) = l^{x) + hfiuix) + (7 Ata (x) + h m {x))C,u 

+ {lau{x) + ha»(x))C,» + MX) + h a p(x))^ a ^ , v (12) 

Here derivatives of £ are with respect to x not x. The generalized gauge transformation 
is defined as one in which the quantity h(x) is substituted for h(x) into the expression of 
interest. If h, £ and their derivatives are small enough, then to leading order this gauge 
transformation reduces to the usual transformation (9). 

It is proven in Ref. |J that, when the vacuum field equations are satisfied, AG^ is 
invariant under this generalized gauge transformation. It is also argued in that paper that 
the quantity 

AiG^ + ... + A n G p 

is gauge invariant to nth order. This implies that Eqs. (11a) and (lib) are invariant to 
nth order under generalized gauge transformations. Thus there is no problem with gauge 
invariance to any order so long as generalized gauge transformations are used. 

In the same paper it was shown that there is a large amount of freedom available in 
choosing the form of the effective stress-energy tensor for gravitational waves. This freedom 
is related to the freedom one has in choosing in the split between the background and the 
perturbed geometry. Thus, as noted above, the BH ansatz 

< AiG^ >= (13) 

is simultaneously a choice of the form of the stress-energy tensor and a definition of the 
separation between the background and perturbed geometry that is to be used. From the 
above discussion it is clear that the condition (13) is, in general, only gauge invariant to first 
order. The lack of exact gauge invariance in this condition is a reflection of the fact that the 
splitting between the background and the perturbed metric is inherently gauge dependent. 

As mentioned above, an explicit calculation which we have made using the gauge trans- 
formation (17) below shows that, even after time averaging, < A2G7M, > is not separately 
gauge invariant. The problem is simply that the condition (13) was used along with the 
expansion (5) in the derivation of Eqs. (6a) and (6b). However, this condition can only be 
imposed in a particular gauge since it is not exactly gauge invariant. Since BH solved Eqs. 
(6a) and (6b) and since we also use those equations as the starting point of the derivations 
in Sec. IV, it is necessary to show that these equations are consistent with the second order 
version of Eqs. (11a) and (lib). 

To begin note that to actually solve Eqs. (11a) and (lib) in a particular gauge one can 
follow Isaacson [|Hj] and make the following expansion for h 

h = h?> + h® + ... , (14) 

where is defined such that AiG(7, h^) is of the same order as A ra G(7, h^). Then to 
second order Eqs. (11a) and (lib) can be written 
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A 1 G( 7 ,/i {1) ) = < A 1 G( 7 ,/i {1) ) > , (15a) 
AiGfr, /i (2) ) + A 2 G( 7 , = < A 1 G( 7 , h^) + A 2 G( 7 , > (15b) 

G( 7 ) = -< A X G( 7 , + A a G( 7 , /i {2) ) + A 2 G( 7 , h®) > (15c) 

Condition (13) can be imposed by requiring that < AiG( 7 , h {n) ) >= for all n. Then 
to lowest order Eqs. (15a)- (15c) are equivalent to Eqs. (6a) and (6b) with the substitution 

It is important to emphasize that it is the original equations (11a) and (lib) that are 
gauge invariant to nth order and that Eqs. (15a) - (15c) and thus Eqs. (6a) and (6b) are in 
general only correct in a particular gauge. The procedure being followed here is the usual 
one for finding solutions to Einstein's equations. First the equations are written down in a 
particular coordinate system and then they are solved (here in an approximate manner) in 
that coordinate system. 

B. A valid expansion of the Einstein tensor 

Since the gravitational waves making up the BH geon are high frequency, large angular 
momentum waves, it is necessary to find an expansion of the Einstein tensor that is appropri- 
ate for these waves. The key point for a thin-shell geon is, that while the waves' amplitudes 
are much smaller than unity, their derivatives are large compared to unity. Further, while 
the background metric is of order unity, some, but not all, of its derivatives are much larger 
than unity in the active region (although they are of order unity or smaller well outside of 
the active region). Thus, whether a term is of leading order will be determined not only by 
the power of h but also by the power of the frequency u or of the harmonic order I* that it 
contains. 

What must be done is to find a self-consistent expansion of the Einstein tensor with the 
above constraints on the background metric and the perturbations. Since some derivatives 
of the background metric are large in the active region but of order unity outside of it, this 
is a complicated task. It clearly depends on the solutions to both the wave equation and 
backreaction equation which in turn depend on the expansion used. 

Fortunately, as is shown by direct calculation in Sec. IV, to second order the method 
used by BH and discussed in Sec. II works. However for the reasons discussed above, it 
would be significantly more difficult to compute the third and higher order terms in the 
correct expansion of AG. 

IV. THE WAVE AND BACKREACTION EQUATIONS 

Derivations of the wave and backreaction equations have been attempted at the time of 
the BH paper || and subsequently by us and independently by ||. Our version provides 
a correct, independent treatment of these important equations. We begin with the wave 
equation. 
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A. The wave equation 



The geon is composed of a large number of gravitational waves and the equations for 



each are the same. BH used the equations AiR 



)1V 



0, with the Ricci tensor. 



BH 



showed that to leading order these equations are equivalent to (6a). We first expand AR^ V 
in powers of and its derivatives, so the wave equation is the high frequency and large 
angular momentum limit of 



i*[h 



fiv;af3 



+ h, 



h 



- h 







(16) 



Here indices as usual are raised and lowered by the background metric jap. 

Following BH we analyze the angular behavior of solutions to these equations for static 
spherically symmetric spacetimes by the method of Regge and Wheeler [|ll[] , and limit our- 
selves to azimuthally symmetric, odd parity functions. In the gauge of Regge and Wheeler 
the perturbation h is obtained from the general h via a gauge transformation of the form 



f = 



-Vjjt 



d 



(17) 



(r 2 sin ( 



l , and 



Here A(r) is some function of r, is an antisymmetric tensor with e 32 
Yi°(9) is a spherical harmonic. In their gauge the "odd" solutions are then written in the 
form 

dY t \e) 



h 



»>l r J 



-iujt 



sin 9 



do 



c.c. 



(18) 



When these equations are substituted into Eq. (16), and when the background metric is 
of the form (7), there result the following three equations for h and hi'. 



lUJt 



' dh _ 2ho\ h 
dr r I 1 



*2 



2 -v , 

uj e + 



d 2 h c 
dr 2 



+ iuj 



dh\ 
dr 



hi 



-(X r + v r 



iu;h e 
1 



r 
1, 

-\v r 
2 V 



\ T — v r — — 

r 







(A r + v r 



dh, 



o 



dr 



A;*2 



2v r 
r 



hn = 



(19a) 
(19b) 
(19c) 



Here we have put I* 2 = 1(1 + 1), and denoted radial derivatives of the background metric 
by subscripts. From Eq. (19b) it is clear that if we take hi to be real then ho is imaginary. 
Taking note of the complex conjugate term in Eq. (18) it is then seen that the part of the 
disturbance associated with ho has a time dependence of the form sin(a;£) while the part of 
the disturbance associated with hi has a time dependence of the form cos(u;£). 

By using Eq. (19b) the function h can be eliminated from Eqs. (19a) and (19c). Then 
the change of variables 

e^ x ~ u) dr , 



dr* 
Q 



hi 



e 2V ' — 



(20a) 
(20b) 



9 



results in the two equations^ 



d 2 Q 
dr* 2 
d 3 Q 
dr* 3 



+ 



+ 



e 2 



2r 



—A r u 3 

— v r * + e 2 + 5 - 
r 



dr 



*2 



+ e 2 



or - e" 



a_i_^ 3 z/*i * 



2r 







7*2 



— e 2 





3 A r * 
2r 



(21a) 



__2._|_i: 3 ^, 

+ e 2 + 2 



2r 



dr* 



+ 



e 2 



7*2 

2 — r + e 



-A+J2 



2 za. 



A _,_ £ A 

_l_ e 2 +2 _ 



4r 



-A I i> Ay* *</•■* 

-e 2 ^2 



2r 



e 2 + 2 



4r 



— A | i/ J/« * j- * — A ■ v 3 Co' 

+ e 2 2 he 2 2 

2 r r 



g = o . 



(21b) 



To ascertain the correct limiting form of these equations we take a hint from the EM geon's 
background and put 
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I* 



dr 



dX 



« 0(1*) 



dr* 



(22) 



(In the following subsections we show that this indeed allows self-consistent solutions.) To 
leading order (/* 2 ), Eqs. (21a) and (21b) are 



d 2 Q 
dr* 2 

d 3 Q 

dr* 3 



+ 



e u l 



.21 
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+ 



u 2 - 



e u l* 



dQ 

dr* 



Q = 
= 



(23a) 
(23b) 



Note that these equations are consistent to leading order in /*. Note also that it is not 
being assumed here that dX/dr* and d 2 vjdr* 2 are of order unity or less, simply that they 
are not of order I*. It turns out that for the geon solution both of these quantities are of 
order unity outside of the active region, but they are of order l* 2 ^ 3 inside the active region. 
It also turns out that dQ/dr* is of order /* Q outside the active region, but is of order l* 2 ^ 3 Q 
inside the active region. Thus Eqs. (23a) and (23b) describe the leading order behavior of 
the gravitational waves, both inside and outside of the active region. 



6 By taking the derivative of (21a) with respect to r* it is easy to show that the two equations 
are not consistent unless the background metric 7^ is an exact solution to the vacuum Einstein 
equations. This is because for a general background the lack of exact gauge invariance does not 
allow the simplifying gauge transformation (17), as explained above. However, the approximate 
gauge invariance of the wave equation results in the equations being consistent to leading order in 
/* if the background geometry satisfies the backreaction equations to leading order. 
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B. The backreaction equations 



The next step in deriving the equations that result in a gravitational geon solution is 
to consider the effective stress-energy tensor, Eq. (4), for the gravitational waves. As with 
the wave equation, the easiest way to derive the stress-energy tensor to leading order is to 
consider the linear and quadratic terms in an expansion of the Einstein tensor in powers 
of hfj, u . The large I* limit of these terms is to be taken in the same way as it was for the 
wave equation in the previous subsection. The result is averaged over times which are long 
compared with the gravitational wave period. 

The background metric for the geon solution is static and spherically symmetric, so from 
the backreaction equation (6b) it is clear that the stress-energy tensor should have these 
same symmetries. Time averaging automatically makes it static over the geon lifetime. 
Spherical symmetry is a more difficult problem. Arranging the waves so that the stress- 
energy tensor is spherically symmetric has been discussed in detail by Wheeler in the 
case of the electromagnetic geon. The arrangements are identical for the gravitational geon. 
They are briefly summarized here. 

First, as can be seen from Eq. (18), to have spherical symmetry in the large I* limit it 
is necessary to have many gravitational waves with different angular orientations. Thus we 
put 

V = E(V)< ( 24 ) 

i 

with each (h^i a solution to Eq. (6a) with a different value of lu, a different phase, and a 
different polar axis. If the phases are random and the values of lu, while different, are all 
approximately equal and large, then after the time averaging is done, to leading order in u 
the cross terms in < A^^ > vanish and 

< A 2 G flu > « J2 < ( A 2^) l > . (25) 

i 

It is useful to denote {A 2 G ilv )i as the wave which has a pole at 6 = 0. Wheeler has 
shown that if the distribution of the polar axes of the waves is uniform, then when the waves 
are all added together one has 

< A 2 Gt > = \[ < {^G t t ) i > sin 0d6 

< A 2 G r r > = - r < (A 2 G r r )j > sin 9 dd (26) 

2 J o 

< A 2 Gg e > = \[ {< (A 2 G e e ) i > + < (A 2 G/) i >) sinfldfl . 

Finally there is the question of conservation. From the discussion of Sec. III. A it is clear 
that < 7] t „ > should be automatically conserved to leading order. Since the geon is here 
treated to leading order, that is in the large I* limit, one expects conservation to hold. Direct 
calculation shows that the stress-energy tensor is conserved to leading order in this limit. 

In the large /* limit the backreaction equations coming from the (t,t), (r, r) and {0,0) 
components of (6b) in an orthonormal frame are 
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— A V Ay* 

e 2 2 

r 



-Z* 



e 2 2 z/_ 



8 (2/ + 1) cu 2 r 4 
8 (2/ + 1) cj 2 r 4 



v, ,2 „2 ,n2 



-r 2 Q 2 + e-W 2 Q 2 



cZr* 



-e 2 2 QrQ- he 2 2 _ Q _ 1_ e v r *\ 

dr* I* dr* \ dr* , 



-(*+") / v r * — X T 
e 2 



2r 



+ e" 



V 



(21 + 1) w 2 r 4 



2/* 2 g 2 + e~ i/ 7^r"g i 
2 



. .2 2 n2 



+ e 



-2i/ 



4cu 2 r 4 fdQ 
dr* 



In arriving at these equations the following identities were used 



1 r l 



2J-i 



[Pi{x)fdx 



1 r 1 



(xPi(x)P/(x)dx 



1 r 1 



(P/(x)) 2 dx 



1 r 1 



(l-x 2 )(P/(x)) 2 dx 



21 + 1 
l 

21 + 1 
2 

l(l + i) 

(21 + 1) 



(28) 



Here Pi(x) is a Legendre Polynomial. 

The combined set of Eqs. (23a), (23b), (27a), (27b), and (27c) completely specify the 
background geometry and the gravitational waves in the high frequency and angular mo- 
mentum limit. In the next subsection Wheeler's expansion is used to show that in and near 
the active region it is possible to arrive at a set of equations which are independent of both 
u and I*. 



(27a) 



(27b) 



(27c) 



C. Wheeler's expansion for the active region 

The active region is defined as the region in which the gravitational waves undergo radial 
oscillations^ rather than radial damping, and where their amplitudes are significant. As can 
be seen from Eq. (23a), the waves will radially oscillate whenever to 2 > e u l* 2 /r 2 . This 
happens at large r, but in order to construct a geon that lasts for a long time, we must 
use waves whose amplitudes are completely negligible there. Radial oscillation also happens 
near r = a, where from Eqs. (8) and (23a) one finds^ 



'By radial oscillations we mean here that the waves oscillate as a function of the radius, not that 
they oscillate in the radial direction. 

8 Until the equations are solved, the location of the active region is unknown. Wheeler [B] initially 
expanded about p = I* . After solving the equations he found the active region was near p = Z*/3. 
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a ~ — 

UJ 



(29) 



By expanding all of the relevant quantities in powers of Z* 1/3 , Wheeler was able to arrive 
at a set of equations which describe the electromagnetic geon in and near the active region 
and which do not explicitly depend on the value of I*. A similar derivation is presented here. 
Before making the expansions it is useful first to change to new variables, 



r = 
-a = 

e" = 

Q = 



P_ 

UJ 

1 - 



2L(p) 



[8(2/ + 1)] 1 / 2 



I* 



-f(p) 



To remove the dependence on I* a new radial variable x is defined such that 

x=(p* - l*)l*~ 1/3 
dx = l*~ 1/3 dp* . 

Then the following expansions are made 

p =r + r 1/3 r (x) + ... , 

L= l*X (x) + Z* 2/3 Ai(a;) + l* 1/3 \ 2 (x) + ... , 



2/3 



/= r 1/3 0(x) + 0i(x) + /*" 1/3 2 (x) + 



i - 



l*- 2/3 j(x)k(x) + ... 



(30a) 
(30b) 

(30c) 
(30d) 



(31a) 
(31b) 



(32a) 
(32b) 

(32c) 
(32d) 
(32e) 



It is next useful to derive several intermediate identities. Using Eqs. (32a)-(32c) in Eq. (32e) 
one finds 



Ai = qik 3 

Combining Eqs. (20a), (30a), (31b), (32a), and (33a) one finds 



(33a) 
(33b) 



He used a scale invariance present in the original equations for the EM geon (and also present in 
Eqs. (23a) and (23b) and (27a)-(27c)) to rescale the solutions accordingly. To facilitate comparison 
with the equations derived by Wheeler we also expand about p = I*. 
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(33c) 



Finally Eqs. (23a), (27a), (27b), (30a)-(30d), (31a), (31b), (32a)-(32e), (33a)-(33c) can be 
combined to obtain, to leading order in I*, the equations 



dj__ 3 + ]_ 

dx k 2 




+ jk<j) = , 








(34a) 
(34b) 

(34c) 



These equations are exactly the same equations that Wheeler found for the electro- 
magnetic geon. Thus we see that, to leading order, the geometry of the gravitational geon 
is exactly the same as that for the electromagnetic geon. This was already found to be the 
case by BH for the spacetime outside of the active region. It has now been shown to be true 
for the geometry inside of the active region. 

This result was hinted at when BH noted that in the large frequency and angular mo- 
mentum limit the equation of motion for the gravitational waves is identical to that for 
electromagnetic waves. It is also consistent with the result of Isaacson that when the 
effective stress-energy tensor for high frequency, large momentum gravitational waves is av- 
eraged over a small region of spacetime, it is of the same form as that for electromagnetic 
waves. 

For a thin shell geon it is necessary that the solution to the wave equation (34a) fall 
off rapidly on both sides of the active region. Equations (34) are therefore similar to an 
eigenvalue problem in which the amplitude of <fr ( as measured, for example, by its peak 
value) plays the role of the eigenvalue: through equations (34b, c) the amplitude determines 
the strength of the "potential" jk, which must be just right to have as a zero energy 
eigenfunction. For the case of the lowest eigenfunction these equations have been discussed 
numerically by Wheeler ||, and more recently by Cooperstock et al |J. The geon metric 
outside the active region, Eq. (8), was first derived from the numerical solution.^ Work is 



9 Cooperstock et. al. || found a different set of equations. Part of the difference is related to the 
fact that they used a different normalization for the variable / in Eq. (30d). This does not lead to 
qualitatively different solutions. However, sign errors in their equations do result in qualitatively 
different solutions. The right hand sides of their Eqs. (4.22) and (4.23) have the wrong signs given 
the sign conventions used in their paper. These sign errors result in overall sign errors on the right 
hand sides of their Eqs. (4.30), (4.31), and (4.46). They also result in a sign error in the last term 
on the right hand side of Eq. (4.47). 

10 Wheeler found e u ~ 0.11, and Ernst [ffl] showed e u = 1/9 exactly. The method of BH yields 
the same exact result for all fields that form geons and whose (effective) stress-energy tensor has 
vanishing trace. Our development in the Appendix yields yet another way to derive the exact 
result. 
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currently in progress to investigate the properties of the other solutions. The results will be 
presented elsewhere. 

Although the numerical results constitute strong evidence for the existence of geon solu- 
tions, without an existence proof of exact solutions to the basic equations (34) the viability 
of any type of geon is open to doubt. || We provide the existence proof in the Appendix. 

A few words are in order about extensions to higher orders in the expansion. In the 
present approximation the waves are a kind of "null fluid" that can be isotropically dis- 
tributed in the active region. If the expansion of the Einstein tensor is carried out to higher 
order in I* the "graininess" of the stress-energy due to non-negligible cross terms of finite 
wavelength in the stress-energy tensor becomes apparent. This invalidates Eq. (25) and 
makes the stress-energy tensor much more difficult to compute. The background geometry 
could however remain spherically symmetric by a different choice of splitting between waves 
and background. 

V. CONCLUSIONS 

We have obtained a correct self-consistent set of equations for the gravitational geon, 
which describe the gravitational waves and the background geometry. These equations are 
accurate in the high frequency, large angular momentum limit. In and near the active region 
they have been shown to be the same set of equations as those found by Wheeler for the 
electromagnetic geon. Thus, to leading order, the geometry both inside and outside of the 
active region of the gravitational geon is identical to that of the electromagnetic geon. 
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APPENDIX: DISCUSSION OF THE GEON EQUATION 

Our geons are taken to be governed by Eqs. (34). These equations focus on the active 
region, and neglect the other oscillating region (beyond r ~ 2.67M, where u 2 > e v l* 2 /r 2 , 
according to Eq. (8)). We can therefore solve Eq. (34a) as a true bound state problem in 
the effective potential jk, rather than a virtual bound state. The appropriate boundary 
condition on (j){x) is therefore 

0(-oo) = 0(+oo) = , (Al) 

From Eq. (34b) it is seen that we can follow Wheeler || and choose the boundary condition 
on k to be 

fc(-oo) = 1. (A2) 
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(Since e u \ r= o = /c 2 | a;= _ 00 this contradicts Eq. (8), and is due to the use of Wheeler's expansion 
before rescaling; see footnote 10. Nevertheless, the ratio e u \ r= o/ e u \ r=QO = k 2 \ x= ^ O0 /k 2 \ x=oo 
will be correctly given.) A third condition follows from Eq. (34c), 

j(-oo) = -oo. (A3) 

We cannot expect to solve Eqs. (34) with any more conditions than these three; the integra- 
tion of the equation itself will tell us whether solutions of the geon type are possible, and 
what their remaining properties are. 

We reduce the system (34) to a single equation by considering k as the independent 
variable and using Eq. (34b) to convert derivatives, d/dx = —<f) 2 d/dk. We will denote 
d/dk by a prime ('). The unknown function j(k) can then be eliminated between the two 
remaining equations. The resulting third order equation is conveniently written in terms of 
a new function 

H = F 2 F" — 2 — Qk 2 where F = 2 , (A4) 

with the result 



— H + kH' — 0. (A5) 

This can be integrated to yield H = Ak, with A a constant of integration. Thus the equation 
to be solved is 

F 2 F" = Ak + 2 + Qk 2 (A6) 

To evaluate the constant A we use the boundary conditions (Al) and (A2) at x — — oo. 
Because F 2 F" = d 2 F '/ ' dx 2 - F~ 1 (dF '/ 'dxf , Eq. (A6) must yield zero at k = 1, hence A = -8. 
Next we use (Al) at x — +oo to conclude that Eq. (A6) also vanishes for fc(oo), so that 

fc(oo) = i (A7) 

In view of Eqs. (30c), (32c) and the rescaling this yields the geon metric (8) outside the 
active region (with M not yet determined in terms of 0). It is appropriate to rewrite Eq. 
(A6) in terms of a shifted variable, u — k — | that exhibits the symmetry of this equation, 

F 2 F" = Qu 2 -\ (A8) 

The boundary conditions (Al) and (34a) now show that 0, F and jk are even functions of 
u, and by Eq. (34b) we can choose x to be an odd function of u. Thus our four boundary 
conditions can be replaced by the more convenient form 

F' = at u = (A9a) 

F^O when u->l/3 (A9b) 

x = o at u = (A9c) 

2jF = 1 at u = (A9d) 
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Of course only Eqs. (A9a) and (A9b) are needed to solve Eq. (A8). 

To show existence of solutions it is enough to concentrate on the interval J:0<m<1/3. 
Starting with some positive initial value F(0) and condition (A9a), we can always integrate 
Eq. (A8) to larger and larger u as long as F(u) remains bounded away from zero. If F(0) 
is sufficiently large, the solution stays positive to u = 1/3. For example, the estimate 



Qu 



1 /3 

| > -§ shows that F(0) > § (^|) « 0.392 is sufficient. We consider 
all solutions that are positive in the interval /, and call them "solutions in J" for brevity. 
Because the right side of Eq. (A8) is negative in /, these solutions are decreasing functions 
in J. Also, if F\ and F 2 are two solutions in / with F 2 (0) > Fi(0), then the difference F 2 — Fi 
is a finite, increasing function in / because it satisfies 



(F 2 - Fx) 




> 0. 



Therefore solution curves do not cross in /, and the "final value" F(l/3) specifies a unique 
solution in I (as does the initial value, F(0)). At each uel the solutions in I depend 
continuously (and monotonically) on initial and on final values. Now consider the greatest 
lower bound (g.l.b.) of the final values. If this were positive it could be lowered, for example 
by integrating from a smaller initial value, for the integration will run to u = 1/3 unless 
F(u) approaches zero at some u < 1/3; but this cannot happen for finite final values because 
F(u) is a decreasing function. Hence the g.l.b. of final values is zero. The corresponding 
limit of solutions F(u) must therefore be a solution that is positive inO < u < 1/3 and 
approaches zero at u — 1/3. This is the desired "eigenf unction" that satisfies the boundary 
condition (A9b).Q 

It remains to be verified that the solution satisfying the boundary conditions in u also 
satisfies the boundary condition in x. This follows from Eq. (34b), du/dx = —F(u). The 
asymptotic form of F(u) near u = ±1/3, 



F(u) - (12) 



1/3 



U 



111 



\u 



nl/3 



can be integrated to yield 



x — > =F 



(96)V3 



u 



2/3 



+ const. 



Thus u = =Fl/3 does correspond to x — ±oo. (By deriving integral relationships from 
the differential equations the same conclusion, as well as the fulfillment of the boundary 
conditions by all the unknown functions, can be established without using asymptotic forms. 
We refrain from displaying these somewhat involved relationships.) 



11 The corresponding "eigenvalue", the g.l.b. of initial values F(0) 
numerically to be 0.3556.... 



! (0), can easily be found 
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